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SUMMARY 




n> 


b*> 


The procedure developed in an earlier paper for the 
computation of coordinates for Vinti's accurate reference 
orbit is slightly extended to include his treatment of the 
zonal harmonic perturbations of the reference orbit. Herein 
is discussed the method for obtaining the position and veloc- 
ity coordinates due to the more important part of the residual 
potential. 
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INTRODUCTION 

Vinti (Ref. 1) has found a gravitational potential for an axially sym- 
metric planet, given by the expression 


V(A V) 


-HP 


P 2 + C 2 7? 2 


which accounts for more than 99*5 percent of the deviation of the Earth's 
potential from spherical symmetry. A method for correcting for the 
effects of any of the zonal harmonics of the residual potential has also 
been devised (Ref. 2). This includes the residual fourth harmonic with 
coefficient J 4 + J 2 2 , and the third harmonic with coefficient J 3 . The 
latter gives rise only to short periodic and long periodic effects without 
singularities. Because of the work of Kaula 1962, and King-Hele, Cook 
and Rees 1963, giving values of J 4 + J 2 2 for the earth ranging from -(0.9) 
10~ 6 to (0.4)10‘ 6 , it has been decided to accept the value J 4 + J 2 = 0 
given by Vinti for purposes of orbit computation. 


1 



ALGORITHM 


1, Sx, Sy, and Sz 

From Ref. 1, we have 

l 3 
n = fj 2 a 2 

(V - 4) = v 0 

g 0 = - M s +2jT ( V 1 - V 2 ) t 

g = - M 

E = M S + E„ + Ej + E 2 , and 
v = M s +v 0 +v x +v 2 . 

From Ref. 2, we then compute; 

the third harmonic short periodic quantities: 

L3 , G 3 , 'tj > g 3 and h 3 , 

the third harmonic long periodic quantities: 

<Xr 'V ^ 

G 3 ’ V S3’ 31x1 h 3 ‘ 

With cr = J 4 + J 2 2 =0, then the variations in the Delaunay variables are 
given by: 

sl=j 3 l 3> sh = o 

J 3 ~ 

W = JjG 3 t -i G 3 

J 2 

u + r \ 

J 2 

Sg = JjS3 t ^ g 3 
J 2 


2 


and 


Sh = J 3 h 3 + r *v 

j 2 


The variations of the elements a, e, and rj 0 are then. 


3a = — SL 
an 


Se = — SL - (ae)' 1 1 SG, and 

fj.e \P-J 


1 - vl 

8 Vo = 2 SG 

0 P V 0 P 


The variations in the uniformising variables E, v, and xp are; 


SE = (p/a) (h'i + sin E Se) 

Sv = (1 - e 2 ) 1/2 (p/a) [SE + (1 - e 2 )' 1 sin E Se] , 

and 

Si// = Sv + Sg 

The variations in the spheroidal coordinates p, 17, and </> are; 

S p ~ (p/a) Sa - a cos E Se + ae sin E SE 
Sp = ( V/V Q ) 8rj 0 + rj 0 cos xp S <//, 
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and 


S0 = Sh + (l - rj 2 sin 2 <//) " 1 (l - ?] 2 ) 1 / 2 


Then 


and 


Also; 


S0 - (MP)' 172 sin 20 SG 


^total ^reference ^ 

^total — ^reference ^ 

^total — ^reference ^ 


‘total = /Kotal+ c2 ) COS <*>totaI 


Vtotal = /(^otal + C2 ) l 1 * Ctotal) sin ^totaX 


and 


Z total ^total ^total 


2. Sx, Sy, and Sz 


Since, as explained in Ref. 2, secular terms are taken through 
order J 2 = J 3 , and long periodic terms through order J 3 /J 2 = J 2 ; then, 
in the determination of Sx , Sy, and Sz, all terms involving a variation 
multiplied by c 2 are ignored. 
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From the expressions; 


where 


and 


P = -f 


h 2 = 


p 2 + c 2 p 2 
/ 0 2 + c 2 


ae J -2a l (p 2 + Ap + B) sin E 
(p 2 +c 2 ) 


We have, 


p 8 p [(p 2 + c 2 ) - (p 2 + c 2 T 7 2 )] + c 2 ?] 8 t 7 (p 2 + c 2 ) 
oh = t = 0 


Also, 


where 


and 


„ fSa Se . - sir 8a l 
8P P =P p |- + ^ + CotESE+ 2^ 

+ 1 \(2p + A) 8 p + p 8 A + SB] 

2 (p 2 + Ap + B) 

2 pSp \ 

( P 2 + c 2 )J 


Sa 


l 


(Sa + 8 b ) 


Sbj = 0, and consequently, 8 A = 0 . 
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Now, 


and 


Also, 


and 


SD = (ap-c 2 ) (aSp + pSa - 2c 2 r} Q 8 

+ (ap-c 2 p 2 ) (aSp + pSa) + 8c 2 (a 2 p 0 Sp Q + p 2 aSa) 
= 2ap (aSp + pSa) 

Sp = (1 _ e 2 ) Sa - 2ae Se 


SB = 



SD' = 8ac 2 


Sa ( J - r lo) ~ a Vo S V 0 + SD = SD 


Therefore, 


SP = P J — + — + cot E SE + — 
P p I a e 2a 


1 


+ 2 p 


2 (p 2 + Ap + B) (p 2 + c 2 ) 


2(p 2 + Ap + B) 



and then, 



Now 


Sp = 



hSP -2PSh 

7) T) V V 

h 3 

v 
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Here, 


Sh = 

v 


2 p (1 - 7 I 2 ) S P + 2p 2 7]hr] 


2\2 


(i-V) 


and, 


Since, 


Sa i 

8P = P *< — - Tan + - — + — -j— rr— 

77 77 1 P 0 2a i (^2 - T i ) 


217877 


(1 - 77 2 ) 


7 ? 2 2 = 


1 

r 


1/2] 

a 2 - 2a j C 2 J 

1 - 

, 8 a i c2 ( a 2- a 3) 

l 

1 

X 

L + ( a 2 - 2a jC 2 ) 2 J 

J 


then. 


and 


877 = 0, 


"877 8 a j ( 2 1 

SP = P J — - - Tan xpBip + — - + - r-z rr 

77 77 I 7 ?,, "I 1(1 -*? 2 ) ( 1 ?2 ” ^ 2 ) 


Also, 


WJ h * 


where. 




= -^- 1(1 - T7 2 ) plp -( p 2 + C 2 ) 77S77 

\ \ 
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Using the expressions, 


x 


y (p 2 + c 2 ) (1 - T) 2 ) cos 4> 


y = -j (p 2 + c 2 ) (1 - 17 2 ) sin 4 > 

and 


z = pn 


We then compute; 


Sx = — - ■ ■ - ((1 -v 2 ) P*P - (/° 2 + c 2 ) 77877) 

V (P 2 + C 2 ) (1 - 7 } 2 ) 


- -/(p 2 + c 2 ) (1 - 77 2 ) sin <£S<£ 


8y = ■ ■■ n 9 ((1 - i7 2 ) pSp - (p 2 + c 2 ) 17877) 

/ (P 2 + c 2 ) (1 -T7 2 ) 

+ / O 2 + c 2 ) (1 -T7 2 ) cos <£§<£ 


and 


8z = + TjSyO 


Then, from the quantities. 
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and 


x = x 


/ P/0 

VV \ 

y a 3 

U 2 + c2 

1 - V 2 ) 


/ pp 

VV \ 

xa 3 

[p 2 +c 2 

l-v 2 ) 

' H 


2= PV + T ) p . 


s* J-A. 

\P 2 + c 


i-v 2 ) 


+ x 


(p 2 + c 2 ) (p 8 p + pSp) - 2 p 2 p 8 p 


(p 2 + c 2 ) 2 


(1 - V 2 ) (17877+77877) - 2 t) 2 t]St} 1 


(I-77 2 ) 


2\2 


h ^, ( a 3 s y + y Sa 3 ) + 2 y a 3 Sh ^ 


h l 


%' [ PP vr} \ S 

Sy = “5 2 - , -1 Sy 

\p 2 + c 2 1 — 77- 


"(p 2 + c 2 ) (p8p + p8p) - 2p 2 p8p 
l (p 2 + c 2 ) 2 


(1 - 77 2 ) (77877 + 77877) - 2xi 2 7]Bri 
(1 - 7? 2 ) 2 

h <*, ( a 3 Sx + xSa 3> + 2xa 3 »»* 
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and, 


Sz = pBr) + 7}Sp + rjBp + pBr) 


From these, we can write; 


i , = x , + Si 

total reference 


^total ^reference ^ ^ 


and 


i , , = z - + Sz 

total reference 


REMARKS 

This algorithm constitutes an addition to the orbit generator (Ref. 3) 
of the Vinti program. It is brought into use only after the Izsak elements 
have been determined by comparison with observation. 
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